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1 Introduction 



One of the simplest models of magnetism is the Heisenberg model describing the ex- 
change interaction of spins. Despite its simplicity, in general only approximate meth- 
ods are available for its study. An exception to this situation is the one-dimensional 
spin-1/2 case. Very early Bethe [1] constructed the eigenstates for the isotropic spin- 
1/2 Heisenberg chain. Much later also the fully anisotropic spin-1/2 XYZ chain was 
discovered to be integrable [2,3] as it is related to the exactly solvable eight-vertex 
model [4]. In this way the spectrum of the XYZ chain and the ground state corre- 
lation lengths are known, cf. [5,6]. The thermodynamics of this model were studied 
in [7] by an elaborate version of the method used in [8]. Unfortunately, only the 
free energy of the XYZ chain could be studied within this approach, the correlation 
functions at finite temperature remained out of reach. 

In this paper we apply an alternative method to the thermodynamics of inte- 
grable quantum chains giving the free energy and correlation lengths. Our treatment 
will follow somewhat the approach of [9,10] where the Suzuki- Trotter formula was 
employed leading to a mapping of quantum chains at finite temperature to classi- 
cal two-dimensional lattice models. The quantum transfer matrix of the Heisenberg 
model was identified as the diagonal-to-diagonal transfer matrix of the exactly solv- 
able six-vertex model which is identical to the row-to-row transfer matrix of an 
inhomogeneous six-vertex model. The eigenvalues are known in terms of a Bethe 
ansatz [11,12,13], the largest one yielding the free energy and the next-leading ones 
the correlation lengths. The study of the limit of infinite Trotter number poses a 
certain problem. In [9] this limit was taken numerically by extrapolation of the 
eigenvalues, in [13,14] it was done analytically for low temperatures. In [10] the limit 
of infinite Trotter number was taken analytically at the Bethe ansatz level for any 
finite temperature. The derived equations are not of the integral equation type as in 
[7]. However, they are very useful for numerical calculations. Unfortunately, only the 
zero temperature limit of them could be studied analytically [10]. This is the main 
motivation for us to choose another approach to the study of the quantum transfer 
matrix. 
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In section 2 we review the two-dimensional eight-vertex model and its relation 
to the spin-1/2 XYZ chain. The quantum transfer matrix is introduced by a direct 
mapping of the XYZ chain at finite temperature to an inhomogeneous eight-vertex 
model. In section 3 the eigenvalue problem for the quantum transfer matrix is in- 
vestigated. The Bethe ansatz equations are analyzed within the 'non-linear integral 
equation' approach [15,16,17,18]. This allows for a very simple treatment of the 
infinite Trotter number limit. For the partially isotropic XX Z chains, also the ther- 
modynamics in an external magnetic field can be studied. As a result we obtain 
integral equations for the largest eigenvalue and the next-leading eigenvalues of the 
quantum transfer matrix. In section 4 certain limiting cases are investigated analyti- 
cally, in particular the low-temperature limit. The critical XX Z chain in a magnetic 
field is studied in section 5. The low-temperaturc asymptotics of the specific heat 
and correlation lengths are derived and compared to predictions by conformal field 
theory. The non-critical XX Z chain and the related spin-1 purely biquadratic chain 
are investigated analytically as well as numerically in section 6. Several technical 
details and calculations are deferred to three appendices. 



2 The quantum transfer matrix of the XYZ chain 



We first consider the symmetric eight-vertex model on a square lattice with peri- 
odic boundary conditions. Each bond of the lattice is occupied by an arrow which 
is pointing either up/down or right/left, respectively. The only allowed arrow con- 
figurations are such that there is an even number of arrows pointing into (and out 
of) each vertex of the lattice. There are eight-vertex configurations with associated 
Boltzmann weights a, 6, c and d, cf. Fig. 1. This model is exactly solvable for all 
Boltzmann weights [4] which are conveniently parametrized by elliptic functions 



a 



c 



b 
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d = +ipH{iX)H 



H 



-tiX + v) 



(2.1) 



with the elhptic modulus k, the crossing parameter A and the spectral variable v 
as new parameters. The theta functions H, 0, and the quarter-periods K, K' are 
defined in appendix A. The 'antiferroelectric' regime of the eight- vertex model in the 
new parameters is given by the restriction 

< A; < 1, 0<X<K', -X<v <X. (2.2) 

For convenience we set the trivial scale of the problem to the value 



where the function h is defined in appendix A. The Boltzmann weights ( p.l|) possess 
the socalled 'crossing' symmetry. A rotation of the vertices by 90° leads to an inter- 
change of configurations 1, 2 with 3, 4 thus amounting to an interchange of a and b 
or a substitution of v by —v. This relation is depicted graphically in Fig. 2 and will 
play an important role below. 

We next consider the row-to-row transfer matrix T{v) as a function of the spectral 
parameter v, holding k and A fixed. T{v) is a family of commuting matrices [4] with 
two special points v = ±A at which T reduces to left and right shift operators which 
are exponentials of the momentum operator P. At these points the Hamiltonian 
limit of the model can be taken by the logarithmic derivative of the transfer matrix 

HxYZ= const. ±(lnT)'(±A), (2.4) 

leading to the XYZ chain describing the exchange interaction of 1/2 spins 

1 ^ 

HxYZ = -1^ E (Jxcrfaf^, + JyaJaJ^, + Jz^rfaf^,). (2.5) 

^x,Y,z ^YiQ Pauli spin matrices and we have cyclic boundary conditions Ui+i = ai. 
The interaction coefficients are functions of k and A 

1 -|- A; snh ^A 1 — k snh ^A cnh Adnh A 
"^^^ 2snhA ' "^^^ 2snhA ' = 2^^±^' ^^'^^ 
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where snh , cnh , dnh are related to standard elliptic functions, cf. appendix A. In 
place of ( p.4|) we shall utilize more directly the relation 

T(t(A -v)) = e±*^ - + ^(^') , (2.7) 

with H = HxYZ (up to an additive constant) from which we immediately conclude 



N/2 

= ,-PH + Oil/N)_ (2.8) 



N J \ 

This equation relates a product of certain row-to-row transfer matrices to the expo- 
nential exp(— /5i7) where (3 is the inverse of the temperature T and is a large even 
'Trotter' number. The partition function of the Hamiltonian can now be expressed 
in terms of the transfer matrices 



lim Tr 



N/2 

(2.9) 



NJ \ 

So, the partition function of the quantum chain at finite temperature is given by the 
partition function of an inhomogeneous eight- vertex model with alternating rows, cf. 
Fig. 3. For the calculation of this partition function the column-to-column transfer 
matrix is best adapted. This 'vertical' transfer matrix (or quantum transfer matrix) 
is nothing but a row-to-row transfer matrix of an inhomogeneous eight- vertex model 
with alternating spectral parameters ±(A — (3/N). This is easily understood from 
the crossing symmetry, cf. Fig. 2. Thus we obtain the representation 

^ = i'?J.TrT-(A-|.-A + |), (2.10) 

where T(., .) denotes the inhomogeneous transfer matrix with length whose eigen- 
values satisfy Bethe ansatz like equations [4]. The free energy per site is / = 
IrnvL^aoF/L = —1/ l3\im.L^^\\iZ/ L where first the limit N ^ oo has to be taken 
and then L — oo. We are allowed to interchange these limits due to the theorems in 
[19,20]. Standard reasoning then yields 

/=--^Jim_lnA^ax, (2.11) 



where A^ax denotes the largest eigenvalue of the inhomogeneous transfer matrix 
T(., .). All other eigenvalues of T(., .) are separated by a gap which is finite for finite 



temperature even in the limit ^ oo. This is an important difference to and an 
advantage over the transfer matrices on the left-hand-side of (|2.8|) . The next-leading 
eigenvalues give the exponential correlation lengthsQ of the equal time correlators at 
finite temperature 



1 



lim In 



A 



(2.12) 



Lastly we want to comment on the study of thermodynamics of the quantum chain 
in the presence of an external magnetic field h coupling to the spin S = Z)j=i Sj, 
where Sj denotes a certain component of the jth spin for instance Sj. Of course this 
changes ( |2.8|) only trivially 

-| N/2 

(2.13) 



A + - T A-- 



A^ 



^f3hS 



-P{H-hS)+0(l/N) 



On the lattice, the equivalent two-dimensional model is modified in a simple way by 
a horizontal seam. Each vertical bond of this seam carries an individual Boltzmann 
weight e^f^^/"^ if Sj = ±1/2 which indeed describes the action of the operator 



(2.14) 



Consequently, the vertical transfer matrix is modified by an h dependent boundary 
condition. It will turn out that these modifications can still be treated exactly if 
the additional operators acting on the bonds belong to symmetries of the model. 
Therefore, a magnetic field cannot be studied unless limiting cases of the XYZ chain 
are investigated such as the XX Z chain. 



3 Bethe ansatz equations and non-linear integral 
equations 

In [4] the treatment of the eigenvalue problem for inhomogeneous transfer matrices of 
the eight- vertex model was described. We content ourselves with quoting the relevant 
*The derivation is standard involving a spectral representation and the observation that corre- 
lation functions of operators acting on vertical bonds are identical to correlations of operators on 
horizontal bonds (for spectral parameter v — s- ±A). 
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formulae. Each eigenvalue of T(., .) is given by a function A(x) at a; = 0, where A(a; 
satisfies the equation 

1 



K{x)q{x) = 
is a 'known' function with definition 

= <|)i(x)<l>2(x), $1/2(3; 



- iX)q{x + i2X) + iu(^{x + iX)q{x - i2X) 



(3.1) 



ph ( -(x=Fixo, 



Ar/2 



(3.2) 



where h{x) is defined in appendix A, is a constant, s = ±1, and Xq involves the 
temperature via 



xo = A--. 

The function q{x) is 'unknown' with functional form 



N/2 

q{x) = exp(— irx) ]^ h 



(3.3) 



(3.4) 



and zeros yj henceforth called Bethe ansatz numbers which have to be determined 
under the condition that A(x) is analytic for all x in the complex plane. This imposes 
that each zero i/j of the left-hand-side of (|3.1| ) is also a zero of its right-hand-side thus 
yielding the well-known Bethe ansatz equations. 

For the XYZ chain with zero magnetic field we have u = 1 and a sum rule 

N/2 

y. = (1 _ rs + N)K + -(s - 1 + N)K' + 4pK + i2p'K', (3.5) 
i=i 2 

where s = ±1 and r = ±1 are quantum numbers of the system and p, p' some integers 

related to r 

r = ^(s-l + Ar + V). (3.6) 

For the critical and non-critical XXZ models the corresponding limiting cases of 
modulus k —>■ 1 and — > describe the problem. In any case r = 0, but there may be 
less than N/2 Bethe ansatz numbers such that the product on the right-hand-side of 
( p.4| ) runs over less than N/2 factors. Now an external magnetic field h may be 
introduced leading to a non-trivial u namely 



UJ 



^Ph/2 



(3.7) 



A compact way for writing down the Bethe ansatz equations is p{yj) = — 1 where 
p{x) is the ratio of the two summands on the right-hand-side of (|3.1| ) 

' <l>{x + iX)q{x - i2X) ^ ' 

Instead of deahng with coupled non-hnear equations for the Bethe ansatz numbers 
we work with ( |3.1| ) on the basis of functional equations in analogy to the approach 
in [16]. Shortly, the usefulness of the following functions will become obvious 

^i^) = — — TT) 2l(a;) = 1 + a(x), 
j)[x — iX) 

a(x) =p(x + iA), = 1 + a(x), (3.9) 

where a, a etc. denote independent functions which are related by complex conjuga- 
tion actually only if cJ = Xjuj. Quite generally a(x), a(x) are very small for real x in 
the limit of large values of A^. The functions 2t(x), 2t(x) in turn are close to 1. In 
terms of these functions A(a;) can be written in two ways 

ijj q{x — %X) 

A{x + tX) = scu<l>(a; + z2A)|||^^2l(x). (3.10) 

Now we recall that all eigenvalues A(x) are analytic in the whole complex plane. 
For a while we focus our study on the two largest eigenvalues with only real Bethe 
ansatz numbers. Each of these functions A(a;) is non-zero in a strip — 2A < Im x < 2A 
where the logarithm defines an analytic function. The Fourier transform of InA(x) 
can therefore be calculated in two different ways by integrals with integration paths 
along Imx = —A or Imx = A using the two different representations ( 3.101 ) of A 



in terms of the functions q, 2t and 21. Of course both calculations yield the same 
result due to Cauchy's theorem thus imposing a non-trivial relation on the Fourier 
transforms of lng(a;), ln2l(a;) and ln2l(a;). In other words, the Fourier transform of 
Ing(x) is given by the transforms of ln2t(a;) and ln2t(x). Next we observe that the 
transform of Ina(x) is given in terms of the transforms of Ing(x) or, alternatively 
by ln2l(a;), ln2t(x). (In fact this reasoning should be applied to the derivatives of 
all the mentioned functions, because Ing(x) does not admit a Fourier transform in 
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contrast to [Ing(x)]'. For more details particularly concerning the regimes of validity 
of the Fourier transforms see appendix B or the exposition of similar calculations 
in [16,21].) To the final equation for the Fourier transform of In a(x) the inverse 
transform can be applied resulting in the non-linear integral equation 



N 



fcisn ^{x- Xo)sn + Xq] 

ZK ZK 

2K 



In a(x) = — In 

l'2K 

+ / k(x — y)\nQi(y)dy 

J~2K 



2K 



k{x — y — i2\ + ie) ln2t(?/)d?/ 



(3.11) 



where e is an infinitesimally small number. The function sn is the standard elliptic 
sn function (cf. appendix A), here however to a new modulus ki which is defined by 
the requirement that the ratio of the corresponding quarterperiods Ki, K[ is 

t = K- 

The kernel of the integral equation k{x) (not to be confused with the modulus) is 
defined by the series 



k{x) 



1 

AK 



1 ^ sh {k{K' -2X))ch (ikx) 

2 ^ ^ 



^j^^ ch(A;A)sh(A;(K'- A)) _ 
where we introduced the short-hand notations 



TC 



ch (x) = cosh —r^x, 
2K 



IT 



sh (x) = sinh x. 

2K 



The remaining integration constant D has the value 



D = \nruj — 2rA, 



(3.13) 



(3.14) 



(3.15) 



where r = ±1. For a{x) we obtain an equation very similar to ( p.ll|) . Both equations 
close since a and o are related to 21 and 21 by ( p. 91) . The integral equation (|3.11 



is valid for the XYZ model with repulsive (0 < A < K'/2) and attractive {K'/2 < 
A < K') interaction. In the latter case the integration paths of the first and second 
integral have to be shifted into the upper and lower half plane by ±z(2A — K'). For 
this situation an alternative set of equations is more suitable, cf. appendix C. 

Lastly we want to find an equation for A(x) in terms of 2l(x) and 'Qi{x). For this 
purpose we may use any of the two equations ( |3.10| ) as we can calculate In q{x) from 
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2t(a;) and 2t(x). However, it is more convenient to multiply both equations ( |3.10|) 
with the result 



A(a; - iX)A{x + iX) = <l>(x - i2X)<^{x + i2X)^{x)Qi{x) 



(3.16) 



where the right-hand-side does not involve q{x) anymore. This equation is the so- 
called 'inversion identity' for finite A^. Taking the logarithm, it is solved immediately 
by Fourier transforms 



InA(x) 



Inr + AT £ '^^^,^i":^^l^"!tS^'^°^ch {k{K' - 2A)) 



2K 



k=l 



kch{kX)sh{kK') 



+ I c(x - y) \n[miy)]dy, 

-2K 



where c{x) is defined by 



^^^^ = AK 



1 ^ch {ikx) 
2^^^MkX)\' 



(3.17) 



(3.18) 



This function is related to the sn function with modulus ki ( |B.8| ) and other elliptic 
functions (p.9| ). 

It remains to study the limit of infinite Trotter number A^ — oo in ( pj.llp , ( pj.lT) ). 
Using ( |A.5| ) and ( |B.8| ) we first obtain 



In a(x) 



—27Tj3c{x) + In ruj 

+ I k{x — y)\n ^{y)dy — / k{x — y — i2X + it) In 2l(?/)d?/, 



2K 



2K 



(3.19) 



where we have used r = for the two largest eigenvalues. The analogous equation 
for a is obtained from ( 3.19| ) by the replacement of a, 2t by o, 21, and i, Inrcj by — i, 
— Inrcij, respectively. From ( |3.17| ) we find (as we are only interested in A = A(0) we 
set a; = 0) 

, , , Ti(5 ^s\iikX)c\iikiK' -2X)) ^ ^ _^ 

lnA = lnr + ^5I — \Jn.\\l n.T^,\ + / c(x) ln[2l2t(x)]da;. (3.20) 



k=l 



c\i{kX)s\i{kK') 



2K 



Equations ( p.l9| ) and ( p.20|) determine the two largest eigenvalues of the quantum 
transfer matrix corresponding to r = ±1. The leading eigenvalue (r = +1) yields the 
free energy according to (|2.1! 



Pf = Peo 



2K 



2K 



c{x) ln[2l2t(a;)]dx. 



(3.21) 



where the ground state energy cq is given by the series in ( |3.20| ). The next-leading 
eigenvalue (r = —1) determines the length ( p.l2|) of the leading correlation function. 



Numerical solutions to these equations compare well to numerical results in [10]. 
However, the mathematical equations therein are not related in an obvious way to 
our ones. 

All other solutions to ( |3.1|) are obtained by allowing for complex Bethe ansatz 
numbers. For the next-largest eigenvalues of the quantum transfer matrix most of 
the Bethe ansatz numbers are real apart from one number on the horizontal line 
Imx = K' (1-string) or two complex conjugate numbers at an approximate distance 
of 2A (2-string) [6]. The main modification of the properties of the eigenvalue function 
A(x) is the occurrence of two zeros 6i, 62 in the strip — 2A < Imx < 2A. For these and 
also more general situations (|3.19|) and ( p.20| ) are still valid, however with a modified 
integration path C encircling certain singularities of the involved functions [21]. 

For a 1-string consisting of a Bethe ansatz number i/q (Imyo = K') the path C 
follows the real axis with loops in the upper half plane encircling 9i+iX, 62+iX (simple 
zeros of 21) and yo—iX (simple pole/zero of 2t/2t) clockwise. The loops can be removed 
using Cauchy's theorem. Technically this is achieved with the differentiated version 
of ( |3.19| ) (with integration path C). This ensures that the integrands [InSt]', [In 21]' 



possess simple poles at the before mentioned points whereas In 21, In 21 have non- 
isolated singularities. The application of Cauchy's theorem produces contributions 
by the kernel k{x). Integrating again we obtain 

lna(a;) = — 27r/3c(x) + In rtu — — — 

2K 

-K{x - {9, + ^\)) - K{x - {02 + ^\)) + In (^^(^ ' + ^^)) 

h2{x - [yo - tX)) 

+ I k{x — y)\n^{y)dy — / k{x — y — i2X + ie)\nQl{y)dy, 

J-2K J-2K 

(3.22) 

where 

I ^ Mk{K'-2X))M^kx) 
K [X) - ^^x + ^^^^^^ _ , [6.Z6) 

satisfying K'[x) = 27rz k{x), and we have used r = 71 /AK. The function /i2 is defined 
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by 



h2{x) = h (^i 



El 
2K 



(3.24) 



in terms of the function h (cf. Appendix A), now however with an elhptic modulus 
/c2 with quarterperiods 

KL K' - \ 

(3.25) 



1_2 
K2 



K 



The function /i2 enters ( |3.22|) because of an apphcation of ( p.6|) . The integral equa- 
tion ( p. 221) has to be solved under the subsidiary conditions a(6'i,2 + ^A) = a{yo + i\) = 
-1. 

In order to derive the counterpart to ( |3.2CI| ) we have to go back to ( |3.17| ) (with 
integration path C). Taking the derivative, removing the loops of C, collecting the 
contributions due to the residues at 9i + i\, O2 + i\, integrating again we find 



InAi(x) 



In Ab(x) + ln(— r) + In 
+ 



kisn ^{x- ^i)sn - ^2) 



2K 



c{x-y) ln[2t2l(y)]dy, 



(3.26) 



where InAb denotes the 'bulk' term in ( p. 171) . Taking now x = and ^ 00 we get 

,'2K 



In Ai = — /3eo + In r + In 



fci sn t^isn 60 

2K 2K . 



+ 



2K 



c(x) ln[2t2t(a;)]dx. (3.27) 



Quite similarly we proceed with a 2-string consisting of two complex Bethe ansatz 
numbers y± with \m.y± close to ±A. Here C follows the real axis with loops encircling 
6*1 + 2A, 6*2 + i\ (simple zeros of 2t), ?/+ — i\ (simple pole/zero of 21/21) clockwise, and 
y- + i\ (simple zero/pole of 21/21) anticlockwise. The loops around 61 + iA, 6*2 + iX 
and y+ — iX are removed. Integration leads to 

Ina(x) = — 27r/3c(x) + ln(— rcij) 

i^f fa , ■\^\ irf (Q , •\^^ , 1 ^2(3^ - {y+ + iX)) 

-Kix - [Oi + iX)) - Kix - [62 + iX)) + In — -— 

h2[x- [y+ -tX)) 

+ J^k{x- y) ln2l(?/)d?/ - J^k{x - y - i2X + ie) ln^{y)dy, (3.28) 

where C is an integration path along the real axis below y± =1= iX, and we have 
used r = 0. This integral equation has to be solved under the subsidiary conditions 
a(^i,2+^A) = a{y±+iX) = —1. The corresponding eigenvalue of the quantum transfer 
matrix is 



lnA2 = — /3eo + Inr + In 



+ / c(x)ln[2l2l(x)]dx. (3.29) 
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4 Limiting cases 



In this section we study several cases of the XYZ chain which can be treated analyt- 
ically. First we consider the XY chain which is equivalent to free fermions. Second 
we investigate the low- and high-temperature behaviour of the general XYZ chain. 

a) The XY chain 

In this case the crossing parameter is given by the decoupling condition A = K' /2 
where Jy = ( p.6| ). For this value (|3.19| ), ( |3.22| ) and (|3.28|) simplify considerably as 



they are no integral equations anymore {k{x) = const.), but explicit expressions for 
the function a 

In a{x) = -2'n(3c{x) + Inrs. (4.1) 

The two largest eigenvalues of the quantum transfer matrix are then exactly given 
by 

InA = -/3eo + Inr + / c(x) ln[2l2t(x)]dx. (4.2) 

J-2K 

Observing that the elementary excitation energy e(x) and momentum p{x) (as func- 
tions of the rapidity x) are related to c{x) by 

e(x) = ^p{x) = 2ttc{x), (4.3) 
da; 



we find the familiar expression 

1 

I-tt/2 

The correlation length of the (a^a^) correlation function is 



(3f = (3eo-- dpln[l + e'''^^'')]. (4.4) 

TT J-7r/2 



1 1 r/"^ 

dp In 



TT J-7r/2 



1 _ p-/3€(p) 



(4.5) 



Of course, also the results for 1- and 2-string excitations ( |3.26| ) and ( p. 291 ) are exact 



with (|4.1| ) (under the subsidiary condition a{6i^2 + i^) = —1)- These eigenvalues yield 
the correlation lengths of the (a^a^) and (a^ay) functions. 

b) Low and high-temperature asymptotics 
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Next we study the general XYZ chain at low temperatures. The integral equa- 



tion ( p. 191) does not admit an analytic solution, however an iterative procedure is 
conceivable. For such an approch ([4.1| ) takes into account just the lowest order in 
temperature. In this approximation the asymptotics of the free energy is given by 



an equation similar to ([4.4] ). Performing a saddle point integration we find explicitly 

/5/ = f3eo - e-^^« [BiT^^^ + EsT^^/^ ^ o{T^/% (4.6) 

where 

Ki, , , f K l-kf f K \'/' 

^o = 7^{l-h), B,= i , M , B2- " ^ 



2K^ ' \7iKi h J ' " 47rV2 \Kiki{l - k 



(4.7) 

Within the same approximation we obtain the correlation length of the (cx^a^) func- 
tion 

^ = 2e-^^°[5iTi/2 ^ ^^2^3/2 ^ o{T^^% (4.8) 
In addition to the terms in ( [4.6| ) and ( [4.8| ) there will be 0(e~^^'^°) contributions from 



the first iteration of the integral equations. This is not studied here. However, for 
the XYZ chain with attractive interaction {K'/2 < A < K') higher order iterations 
may lead to contributions which become dominant for sufficiently strong attraction. 
For Im a; ~ 2A — i^T' we find 

In a{x) = -27ipc{x) + exp {-27ipc[x - i2{2\ - K')]) + ... . (4.9) 

Within this accuracy we get 

(3f = (3eo — 2 / c(a;) exp (— 27r/3c(x)) — / Cb(x) exp (—2TT/3cb(x)) 

J~2K J-2K 

= r'"^ dpe-^^(p) - — r dpe-'^''^P\ (4.10) 

TT J-n/2 2n J-n 

with 

Cfe(x) = c{x-i{2\-K')) + c{x + i{2\-K')), 

eb{x) = -T-Pbix) = 27rcbix), (4.11) 
ax 

where eb{p) is the energy-momentum dispersion of the lowest bound states [6]. The 
result of ( [I.IU] ) is physically very intuitive. For certain interaction parameters the 
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gap of the bound states may be lower than that of the free states thus becoming the 
thermodynamically relevant one, cf. also [7,10]. 

Without derivation we communicate our result for the low-temperature asymp- 
totics of the correlation length of the (a^a^) and (a^a^) functions 

The second term in the asymptotic expansion of the (a^a^) function decays expo- 
nentially with length ^'^ 



1 , , /27cKV 



Of course, these results are only valid in the repulsive regime < A < K'/2. In 
the attractive regime the lengths are dominated by bound states. For the general 
XYZ chain (eo > 0) the length diverges at T = with an essential singularity 
while all other lengths ^x,y, remain finite. This implies long-range order at zero 
temperature. 

We conclude this section with a comment on the high-temperature asymptotics 
of the free energy. For small values of (3 (and r = -|-1, uj = 1) the function a{x) 
becomes independent of x and ( p.l9[ ) turns into a simple algebraic equation. The 
result is o = a = 1. The integral in (|3.21|) can be done, yielding 



/--Tin 2, (4.14) 



with a high-temperature entropy In 2 as it should be for a model with two states per 
site. 



5 The critical XX Z chain 

Next we study the critical limit of the XYZ chain with modulus — > 1 (entailing 
K oo, K' 7r/2) where (|2.6|) reads 



Jx — — , -Jy — — , 'Jz — — : , 

sm 7 sm 7 sm 7 
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with 7 = 2A G (0,7r). As we are allowed to change any two signs of the interaction 
coefficients due to a unitary transformation [4] we see that the limiting case k ^ 1 
corresponds to the XXZ chain with equal Jx and Jz coefficients and a coefficient 
Jy which is smaller in absolute value. The limiting cases of 7 ^ and vr yield the 
antiferromagnetic and ferromagnetic Heisenberg chain, respectively. For < 7 < 7r/2 
the excitations are due to free states only, for 7r/2 < 7 < vr there are also bound states. 
The dispersion relation of the free states is 



TT 



e(k) = — Isin fcl , 

7 



with an obvious sound velocity of 



vr 
7' 



(5.2) 



(5.3) 



At zero temperature the model is critical, all correlation lengths diverge like ~ 1/T 
which will be shown below. 

The A; ^ 1 limit of (|3l9|) and ( [OOD reads 

7r(3 1 



In a(x) 



7 cosh -X 



+ 



IT 



Incj 



7r — 7 



+ 



and 



where 



k{x — y) ln2t(?/) — k{x — y — i'y + ie) ln2t(y) dy, 
\n[m{x)] 



A = -peo + 



27 



cosh -X 

7 



-dx. 



k(x) — > — 



sinh — 7^ kcos{ka 



00 2 cosh ^/c sinh ^^/c 



-dfc, 



and we have used 



c(x) 



27 cosh -X 

I 7 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



One of the subtleties in the derivation of ( ^.41) is the 'renormalization' of the additive 
constant due to the contribution of the constant term in k{x) ( |3.13| ) in the limit of 
K ^ 00. This can be checked directly by considering the asymptotic behaviour of 
the two sides of (|5.4| ) for x — 00 (Ina(oo) = 21nu; etc.). It appears that the critical 
limit of ( p.l9|) for r = — 1 becomes inconsistent. The missing eigenvalue is described 
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in terms of (pure) hole excitations with N/2 — 1 Bethe ansatz numbers, N being the 
Trotter number. 

In general the excited states are described by equations identical to (|5.4| ), however, 
with more complicated integration paths C Proceeding as in section 3 or simply 
taking the limit of (|3.22|) , (|3.26|) we find for 1-string excitations 

max = 



7 cosh -X 

' 7 



+ m 



In to" 



IT 



7 



sinh^(x- (yo + n/2)) 



+ 



and 



where 



A = -(3eo + In 



K(x) 



k{x — y) In 2l(?/) — k{x — y ~ i'j + ie)\ia 2l(|/) dy, 

ln[2l2t(x)] 



TC TC 

tanh — 6i tanh — 62 
27 27 



1 



27 Jc cosh -X 



dx. 



sinh (f — 7) ^ sin(A;a;) 
2k cosh ^k sinh ^^fc 



dfc. 



(5.8) 
(5.9) 

(5.10) 



For 2-string excitations the equations are identical upon the replacement of yo by 
Excitations corresponding to N/2 — 1 Bethe ansatz numbers are described by 
similar equations where the yo dependent term is replaced by the additive constant 
^^E^7- The equations for the largest eigenvalue for the critical XX Z chain in zero 
magnetic field have also been derived in [18] and [22] using the non-linear integral 
equation approach to the study of Bethe ansatz equations. 

We do not dwell any longer on these integral equations as we aim now at an analy- 
sis of all eigenvalues in the low-temperature limit using an approach with reference to 
[21]. There the finite-size corrections to all eigenvalues of the six- vertex model trans- 
fer matrix were calculated analytically. The integral equations in [21] are identical 
to ( ^.4|) and those for the excitations provided we exchange the 'inhomogeneities' 



TT/3 



7 cosh -X 



N In tanh 



TCX 



(5.11) 



where now denotes the system size of the (homogeneous) six- vertex model in [21]. 
The functions (|5.11|) have the same asymptotic behaviour if we identify 



vl3 = N. 



(5.12) 
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In [21] the 1/A^ corrections to the eigenvalues were found just from the asymptotic 



behaviour of the function on the right-hand-side of (|5.11|) . This means that the 
l/(3{= T) contributions to the eigenvalues of the quantum transfer matrix can be 
obtained from [21] via the correspondence ( |5.12| ). The result is 

In A = -Peo - '^T(x - c/12) + o{T^) + iPo, (5.13) 

with 



and 



x = b H — ; T^Tn +1 — , 5.15 

2 2(l-7/7r) 7r-7 2r' ^ ' 

where the magnetic field entered via ( |3.7| ), 5, m are integers, and Pq = {S — m)-n is 
a lattice momentum. A condition for the validity of ( ^.13|) is the finiteness of a;, i.e. 
h = 0{T). 

Next we discuss the physical consequences of ( |5.13| ). Applying ( ^.111 ) to the 
largest eigenvalue (x = 0) we calculate the free energy density at low temperatures 
and magnetic fields 



/,r..)^e.-i^(r^.^^), (5.16) 



which is consistent with the scaling prediction f{T) = /(O) — (7rc/6f )T^ derived from 
conformal invariance [23]. The central charge c (not to be confused with the related 
quantity (|5.14| )) of the underlying field theory is 1. The low-temperature behaviour 
of the zero-field susceptibility is 

X=7rj^—y (5.17) 
2v[7r — 7j 

The asymptotic behaviour of the correlation functions is determined by the next- 
largest eigenvalues ( |2.12|) with non-zero x 

Cr ~ ~ cos(Poi?)e~^^^^. (5.18) 



The 2-point function for spin components in axial direction is described by 5* = 0, 
m = ±1. For planar spin components we have the selection rule S = 1, m = 0. For 
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non-zero magnetic field the quantity x is complex if the quantum number m 7^ 0. 
The dominant behaviour of physical correlation functions is determined by a pair 
±m, such that 

^ ( h \ _2^Re (a;)Ti? 

Cr ~ cos m K \ e y ^ ' , (5.19) 

Y n — ■y V J 

where we have ignored for the moment the alternating sign for a lattice moment 
Pq = n. The exponential decay of the correlation functions can be obtained from 
conformal invariance [24,25] 

Cr ~ \ e V , (5.20) 



but not the prefactor whose oscillations in general are incommensurate with the 
lattice. Obviously, the scaling dimensions x of the underlying field theory are given 
by X = Re (x) in accordance with finite-size calculations, see for instance [21] and 
references therein. 



6 The non-critical XX Z chain 

Lastly we study the modulus k —>■ limit of the XYZ chain where K —>■ 7r/2, 
K' ^ 00, and fE.6|) reads 



■^^-2ihnrA' '^^-2ihnrA' ^^-~2ihnrA' ^^-^^ 
where A is any positive real number. In this case again two interaction coefficients 
are equal. The third one is larger in absolute value. This limit of the XYZ chain is 
the non-critical antiferromagnetic XX Z chain. The ferromagnetic XX Z chain can 
also be obtained in the limit k ^ 0, however with fixed K' — A, cf. ( |C.1|) . 



The thermodynamics of this model are given by ( p. 191) , (|3.20|) for the antiferro- 
magnetic regime, and by ( P.IO ), (|C.11| ) for the ferromagnetic case where the limit 



k ^ can be performed in a straightforward way. There are two largest eigenvalues 
of the quantum transfer matrix which are exponentially close in the low-temperature 
limit. From these eigenvalues the free energy and the leading correlation length ^ 
can be calculated. Upon approaching zero temperature the length ^ diverges herald- 
ing the onset of long-range order. All other lengths remain finite. In the limit of 
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the one-dimensional classical Ising model (A — > oo) the integral equations turn into 
simple algebraic equations from which the well-known results are recovered easily. 

In Fig. 4a the specific heat and $, are shown for anisotropy cosh A = 3/2, with 
a convenient normalization such that Jx = Jy = 1 and Jz = —3/2. In Fig. 4b 
these curves are given for the same Hamiltonian with ferromagnetic sign. Note the 
divergence of the curves at T = with an essential singularity exp(±A/r) where 
A = I • 0.173... and A = 1 in Fig. 4a and b, respectively. The quantities 0.173... and 
1 are the energy gaps of the excitations for the model with antiferromagnetic and 
ferromagnetic sign. 

The antiferromagnetic XXZ chain of the previous paragraph is related to the 
so called spin-1 biquadratic chain 

N 

^biq = -E(44+ir, (6-2) 

where Sj denotes the vector of spin-1 operators. The relation was first observed in 
[26]. The most systematic way to establish this employs the Temperley-Lieb algebra 
[27] which is satisfied by the biquadratic chain and the previous XXZ chain. The 
ground state energies and gaps coincide. For a calculation of the singlet-singlet 
correlation length at T = see [28] . 

Here we want to describe how the thermodynamics of the biquadratic chain are 
tackled. The problem of the socalled Temperley-Lieb equivalence is the unknown 
degeneracy of each energy level of the biquadratic chain, otherwise its spectrum is 
identical to that of the XXZ chain (for certain boundary conditions). We therefore 
choose a different approach using the classical two-dimensional counterpart to the 
biquadratic chain [28]. This model is an exactly solvable 3-state vertex model. We 
repeat the reasoning of section 2, derive the corresponding quantum transfer matrix 
which is an inhomogeneous 3-state vertex transfer matrix. At this point we make use 
of the Temperley-Lieb correspondence. The two largest eigenvalues of the quantum 
transfer matrix of the biquadratic chain are equal to the two largest eigenvalues of 
the transfer matrix of the XXZ chain 

Abiq(T)=Axxz(T,a;), (6.3) 
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provided that an external magnetic field is applied to the XXZ chain such that 
u! + l/u = 3. (More generally, the quantum transfer matrix eigenvalues for the 
biquadratic chain in a field hb and the XXZ chain in a field h coincide if 1 + 
2 cosh hh/T = 2cosh/i/2T.) The proof of this exact relation will be published else- 
where. In Fig. 5 the specific heat and the leading correlation length is shown. Note 
that the divergence at low temperature is described by the same essential singularity 
as in the a.f. XXZ case. The values of the specific heat in Fig. 5 are generally 
larger than those in Fig. 4a. This is due to the larger high-temperature entropy of 
the biquadratic chain. In this limit the solution to ( |3.19| ) is a = u'^ and a = l/tu^. 
The integral in ( p.21|) can be done, yielding 

f -T \n{u +l/u) = -T In 3, (6.4) 

with a high-temperature entropy In 3 due to the three states per site for the bi- 
quadratic model. The curve for the specific heat shown in Fig. 5 compares quite 
well with the results of [29]. A discrepancy appears at low temperatures since the 
biquadratic model was treated by numerical diagonalization of finite chains in [29]. 

The thermodynamics for the biquadratic chain with opposite sign can also be 
treated, just by introducing a negative temperature for ( |6.2| ). Then equations ( |C.10|) , 
( |C.11|) have to be solved. We do not show numerical results for this case, because the 
computations at low temperature are difficult as it appears that the integration paths 
have to be deformed. However, the low-temperature limit can be studied analytically. 
The result for the two largest eigenvalues is quite interesting 

Amax = -Peo + Ino;, A = -(3eo - Inu, (6.5) 

which implies a residual entropy Sq = Inu = 0.962... and a finite correlation length 
^ = l/(21nco') = 0.519..., with the explicit value to = (3 + \fb)/2. The value of the 
non-zero 5*0 was derived earlier [30]. The result for the leading correlation length is 
novel, its finiteness implies absence of long-range order at T = 0. 



Acknowledgments 

I am grateful to A. Schadschneider and T. Wehner for useful discussions. 



20 



Appendix A 



In this appendix we list the definitions of the Jacobian eUiptic functions and related 
quantities. Elliptic functions are complex, meromorphic and doubly periodic func- 
tions. Here we denote the argument by u. Another parameter which determines the 
periods is the nome q with < g < 1. The quarter-periods K, K', the modulus k 
and the conjugate modulus k' are defined by 

I _ g2n \ 2 



K = 


1 i^/l+g2n-l 


K' = 




k = 




k' = 







l + q 

4 



2n 



(A.l) 

such that k"^ + k'"^ = 1 and q = exp^—nK'/K). 

The Jacobian elliptic functions sn , cn , dn with periods AK, AiK', and the related 
snh , cnh , dnh functions are 

snu = k 2 snhw = —tsmu, 

e{u) 

cnu = — ^. , , cnhw = cnzti, l^-^J 
V y e(xi) 

dn ^ = A; 2 dnhu = dnzii, 

Q[u) 

where the theta functions H{u), Q{u) are defined by 



(A.3) 

These functions satisfy a bewildering variety of identities of which we cite only 
H{u-riK') = i5-V4exp(-||)e(K), 
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H{u) 


^ 1 . TTM 

= 2(74 sm 

^ 2K , 




- n i-2g^ 

n=l ^ 




= H{u + K), 




= e(ii + x). 



e{u + zK') = tq~y^exp(~^)H{u), (A.4) 

V 2K J 

and 

sn (u + iK') = (A.5) 
ksnu 

A more comprehensive summary can be found in [31] and [4]. 

In the main body of the paper a function h{v) is used which is defined by 

h{v) = -t&{0)H{iv)e{iv), (A.6) 

enjoying the properties 

h{v + i2K) = -h{v), h{v + K') = -g"^/^ exp(7rf (A.7) 



Appendix B 



Here we hst the Fourier transforms of some functions appearing in section 3. The 
main strategy of the solution procedure in section 3 is to formulate functional equa- 
tions for certain functions and to solve them by Fourier transforms. All functions 
which we have to deal with are AK periodic. We therefore introduce the Fourier 
transform and its inverse 

oo 

fix) = 

fc=— oo 

r-2K 



1 /■2K 

f, = /(x) e-^^'^Mx. (B.l) 

4a J-2K 



As we have to study the logarithms of periodic functions g{x), it may happen that 
hig{x) is not periodic, but acquires a phase shift lng{x + AK) = lng{x)+ const. In 
that case we apply ( p.l| ) to the derivative of lng{x) or introduce a linear term in 



( [B.l|) . For the functions h and sn we have the explicit formulae 



^ 7 -S- t-r- 



with 



N/2 

7 = 9'/'0(O) 11(1 (B.3) 
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and 



, /•,l/2 \ li -^^ 

^"^^^ / sn X) = - Ing + + _E ^(^^ 



(B.4) 



fc— — oo 



In section 3 we take the Fourier transforms of equations relating different functions 
at different values of the argument. One has to be aware that a Fourier representation 
is convergent only in certain strips of the complex plane. For the largest eigenvalues 
for instance, i.e. real Bethe ansatz numbers, Ing(x) is analytic in each strip n2K' < 
Imx < {n + l)2K' , but singular on the boundaries. We take the strip —2K' < Imx < 
as the 'fundamental' regime where we define the Fourier transform. All expressions 
of Ing(x) with arguments x in other strips are related to the 'fundamental' regime 
by (quasi) periodicity. In this way we find the appropriate versions of ( p.lO ) and the 
appropriate expression for a(x) 



A(x - iX) 



K{x + iX) 



a{x) 



-q'"'' exp f-^(3A - xo)) $i(x - «2A)$2(x + i{2K' - 2A)) 
,(x + ,(A-2V))^ 
q[x — iX) 

swexp f - xo) j + i{2X - 2A:'))'^'2(a; + i2X) 

q{x — iX) 



q{x + i{X-2K')) 
Nn 



2 ~N/2 

ruj q ' exp 



(xo - 3A) 



$i(x) 



q{x — i3X) 



q{x + i{X - 2K'))' 
The kernel (p.l3|) satisfies two important functional equations 



<l'i(x - i2X) ^2ix + i{2K' - 2A)) 

(B.5) 



k{x) + k{x - i2X) 



27T 



In 



h (i^x, k2 



h{i^{x-i2X),k2) 



(B.6) 



and 



k{x) - k{x - i2{K' - A)) = — <^ In 

271 



sn 



(El 

\2K 



X, ki ] sn 



V2a: 



{x-i2{K' -X)),h 



(B.7) 

which may serve for an analytic continuation of fc(x) outside the strip of convergence 
of ( CT) - 
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Lastly we give two expressions of the function c{x) ( p.l8[ ) in terms of standard 
elliptic functions 



cix] 



c{x) 



-I [In sn r I X — I 



:dn 



-xk,] ^-A 



(B.8) 
(B.9) 



where the sn function carries the modulus ki (|3.12| ) and the dn function a new 
modulus k^. 



Appendix C 

In section 3 we have derived a set of integral equations which determine the ther- 
modynamics of the XYZ chain. However, for the model with attractive interaction 
{K'/2 < A < K') the analysis of equations ( |3.19| ), ( |3.20| ) becomes difficult, particu- 



larly in numerical applications as the integration paths have to be shifted into the 
upper or lower half plane. We can avoid this by an alternative approach to the case 
K'/2 <\<K'. From Q we find 



Jx{\) = +JxiK'-X), Jy(A) = -JviK'-X), Jz{\) = +Jz{K'-\). (C.l) 

Recalling that due to a unitary transformation we are free to change any two signs of 
the interaction coefficients [4] we see that the interchange of A with K' — X amounts 
to the replacement H ~H. Instead of treating the thermodynamics of the chain 
with crossing parameter A at temperature T we may study the chain with crossing 
parameter i^' — A at temperature —T. This is actually possible. 

From now on we consider the model for a crossing parameter A with < A < K' /2 
at negative temperatures — T (T > 0). In this case the auxiliary functions a and a in 
are no longer small for large values of N and p. Instead we find that we have 
to take the reciprocals 

a(a;) = p{x — iX), 

a{x) = -^-^—TT. (C.2) 
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In terms of these functions we find two representations of the eigenvalue function A 

q{x — iSX) 



A(x - iX) 
A{x + iX) 



q{x — iX) 



ijj ' q{x + %X) 



(C.3) 



Proceeding similarly to section 3 we derive the integral equation 



In a(x) 



N 

+ 
+D 



In 



sn ¥^{x + i{xo — 2A)) 



2K 



K2 

2K 

sn ^(x — ixo] 



2K 



2K 



k{x — y) hiQl{y)dy — / k{x — y — i2X + ie)ln.^[y)dy 



-2K 



where the sn functions carry the modulus ^2 (|3.25| ) and D is a constant 

N/2 



D 



The kernel k{x) is defined by 



k{x) 



1 
AK 



E 



sh(A;(ir'-2A))ch (ikx) 



\ c\i{k{K' - X))s\i{kX) _ 



(C.4) 



(C.5) 



(C.6) 



and is identical to ( |3.13| ) after a substitution of A by i^' — A. Unfortunately, in this 
case the constant D cannot be expressed in terms of the auxiliary functions. Instead 
of this missing relation we find the subsidiary condition 



2K 



ln2l(x)dx 



AK J-2K AK J-2K 



2K 



ln2l(a;)dx = 2 Ina;, 



(C.7) 



which completes the integral equation. For a a similar integral equation is valid. The 
eigenvalue is given by 

N-K 



InA(x) 



Inr + 



AK 



(xo - A) 



^ ch {k{K' - A)) - ch (^fcx)ch {k{K' - xq)) _ 2A)) 



+ 



k=l 
2K 

-2K 



kch{k{K' - X))sh{kK') 
[c{x -y- i{K' - 2A)) ln2l(i/) + c{x - y + i{K' - 2A)) ln2t(?/)]d?/, 

(C.8) 



where c{x) is defined by 



c(x) 



1 
AK 



+ E 



ch {ikx) 



{c\i{k{K' -X)) 



(C.9) 
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which is identical to ( |3.18| ) after a substitution of Xhj K' — X. In the hmit of infinite 
Trotter number we find 



Ina(x) = -27rf3c{x + i{K' -2X)) + D 

f-2K f-2K _ 

+ / k{x — y) In 2l(y)dy — / k{x — y — i2X + ie) In %{y)(\y, 

J-2K J-2K 

(C.IO) 



and 



In A = In r H 

+ / [c(x + i{K' - 2A)) ln2t(x) + c{x - i{K' - 2X)) ln2l(x)]dx. 

J-2K 

(C.ll) 



These equations determine the free energy density and the leading correlation length 
of the quantum chain in the attractive regime. Here we do not present a study of 
the next-leading eigenvalues due to 1- and 2-strings. 



1 + 2E 



{k{K' -X))c\i {k{K' -2X)) 
ch(A;(ir'- A))sh {kK') 
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Captions 



Fig. 1 The eight arrow configurations allowed at a vertex and the corresponding Boltz- 
mann weights. 

Fig. 2 Graphical representation of the 'crossing' symmetry expressing the invariance 
of the Boltzmann weights of a vertex under a combined rotation by 90° and 
the substitution of v hy —v. 

Fig. 3 Representation of a iV x L square lattice with a distribution of spectral param- 
eters ±vo — ±{X — P/N) homogeneous in each row and alternating in columns. 
N denotes the 'Trotter' number and L the length of the quantum chain. The 
column-to-column transfer matrix is identical to the transfer matrix with row 
inhomogeneity depicted on the right. 

The study of the quantum chain in an external magnetic field leads to the 
introduction of a horizontal seam into the lattice model and modified boundary 
conditions for the column-to-column transfer matrix. 

Fig. 4 (a) Depiction of the specific heat and the leading correlation length for the 
anisotropic Heisenberg chain with Jx = Jy = 1 and Jz = —3/2. (b) Specific 
heat and leading correlation length for the ferromagnetic case Jx = Jy = 1 
and Jz — 3/2. Note that the curves for the correlation length have been scaled 
by factors of 1/100 and 1/1000, respectively. 

Fig. 5 Depiction of the specific heat and the leading correlation length for the spin-1 
biquadratic chain. The values for ^ have been scaled by a factor of 1/100. 
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